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Abstract 

Relay deployment in Orthogonal Frequency Division Multiple Access (OFDMA) based cellular networks 
helps in coverage extension and/or capacity improvement. In OFDMA system, each user requires different 
number of subcarriers to meet its rate requirement. This resource requirement depends on the Signal to 
Interference Ratio (SIR) experienced by a user Traditional methods to compute blocking probabiUty cannot 
be used in relay based cellular OFDMA networks. In this paper, we present an approach to compute the 
blocking probability of such networks. We determine an expression of the probability distribution of the 
user's resource requirement based on its experienced SIR and then classify the users into various classes 
depending upon their subcarrier requirement. We consider the system to be a multidimensional system with 
different classes and evaluate the blocking probability of system using the multi-dimensional Erlang loss 
formulas. 

Index Terms 

Blocking Probability, Inter-CeU Interference, Markov Chain Modeling, Multi-Dimensional System. 



I. Introduction 

The Third Generation Partnership Project Long-Term Evolution Advanced (3GPP-LTE-Advanced) 
proposes different schemes for mobile broadband access in order to meet the throughput and coverage 
requirements of next generation cellular networks [?]. Deployment of Relay Stations (RSs) to increase 
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coverage area and/or improve capacity [?] is one of the proposed techniques in LTE-A. In this paper, 
we analyze the capacity improvement due to RS deployment and analytically determine the blocking 
probability to quantify this improvement. Blocking probability corresponds to the probability that a 
user is denied sevice due to non- availability of sufficient resources in the network. 

Users who experience poor signal strength from Base Station (BS), will require more resources 
to meet their rate requirement and a large amount of resources will be consumed in serving such 
users. This will lead to an increase in the blocking probability. With RS deployed in the network, 
the Signal to Noise Ratio (SNR) experienced by these users improve due to closer proximity of RS 
and now, they will be able to meet their rate requirement with fewer resources. This will reduce the 
blocking probability and improve the system capacity. However, as the radio resources are shared 
between BS and RS, deployment of RSs introduce additional sources of interference. Therefore, it 
is significant to study the impact of interference on the blocking probability. 

Blocking probability has been used as a performance metric in [?], in which transmission scheme 
selection policy (single hop or multi-hop) has been proposed to provide guaranteed target Bit Error 
Rate (BER) and data rate to a mobile user. Another metric to quantify the performance improvement 
in a cellular network is Erlang capacity, which is the traffic load in Erlangs supported by the cell 
while ensuring that blocking probability remains less than a certain value. Despite sufficient work on 
Erlang capacity and blocking probability in cellular networks ([?], [?]), limited literature is available 
on determining the Erlang capacity of relay based cellular Orthogonal Frequency Division Multiple 
Access (OFDMA) networks ([?], [?], [?]). In [?], the uplink Erlang capacity of relay-based OFDMA 
network has been derived considering adaptive modulation and coding supporting both voice and 
data traffic. In [?], the uplink capacity and spectral efficiency of relay -based cellular networks have 
been analyzed. The bandwidth distribution between BS and RSs has been determined to ensure 
blocking probability to be less than a specific threshold. The impact of number of RSs and their 
position on Erlang capacity is investigated by considering Adaptive Modulation and Coding (AMC) 
and Multiple-Input Multiple-Output (MIMO) transmissions. 

In [?] and [?], it is assumed that all users require equal number of resources. However, the impact 
of user location, shadowing and interference from neighboring cells on the resource requirement has 
not been considered. If distinct users of same data rate requirement are present at different locations, 
they may experience different Signal to Interference Ratio (SIR) and hence require different resources 
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in terais of number of subcarriers to satisfy their data rate requirement. In the queuing literature, the 
problem of incoming users requiring different number of resources has been addressed in [?] and 
[?]. In [?], wide band and narrow band traffic is considered, where no queuing is allowed for narrow 
band traffic and a finite length queue is provided for wide band traffic. The blocking probability for 
each traffic is determined using numerical methods. In [?] and [?], the problem of multiple server 
requirement is analyzed and multidimensional Erlang loss formulas have been derived. 

To the best of our knowledge, no literature is available for the computation of blocking probability 
in relay based cellular OFDMA systems, where different users of same rate requirement need different 
subcarriers. In [?], different subcarrier requirement of users has been considered in the blocking prob- 
ability computations. However, SIR experienced by a user and distribution of subcarrier requirement 
were determined using simulations. In [?] (by one of the authors). Cumulative Distribution Function 
(CDF) of interference is computed analytically. However, blocking probability is not determined. 

In this paper, we propose an analytical model to evaluate the performance of a relay based cellular 
OFDMA network (such as an LTE network) in terms of blocking probability. The distinct features 
of our paper are that we consider the impact of user location, shadowing and interference from 
neighboring cells in our analysis for blocking probability, which is not available in the literature. 
Specifically, we determine the SIR experienced by a user and probability distribution of the number 
of subcarriers required. Then, we classify incoming users into different classes based on their 
subcarrier requirement. We consider the network to be a multidimensional system with different 
classes and model the system states by multidimensional Markov chain. In such a system model, the 
computational complexity is more due to the large state space involving the states of both BS and 
RS. To reduce this complexity, we propose an approximation where the state space of BS and RS 
are decoupled. With this simplification, we evaluate the blocking probability of each class in a relay 
based OFDMA network. Admittedly, this is an approximation, however it gives an upper bound on 
the blocking probability. 

The rest of the paper is organized as follows. Section |ll] introduces the system model for the 
downlink of relay based cellular OFDMA network. In Section [nil a model to characterize Inter-Cell 
Interference (ICI) on a Mobile Station (MS) is presented and the CDFs of ICI on BS-MS, BS-RS and 
RS-MS transmission links are derived. In Section |IVl an analytical model is proposed to determine 
the subcarrier requirement and their probability distribution based on ICI experienced. In Section 
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rvl the incoming users are classified into various classes based on their subcarrier requirement. It is 
also shown that complexity is introduced due to the large size of state space when both BS and RS 
are considered. Then, an approximate analytical model is developed by considering the state space 
of BS and RS separately. This model is used to compute the blocking probability for each class 
of user on the three transmission links. Finally, the blocking probability of a relay-based OFDMA 
network is computed using multi-dimensional Erlang loss formulas [?]. In Section |VIl the results 
are discussed and Section IVIII concludes the paper. 

II. System Model 

We consider the downlink transmission scenario in a relay-based cellular OFDMA network as 
shown in Fig{TJ We define the reference cell as a combination of seven sub-cells. The central sub- 
cell Hq consists of a BS centered at (0, 0), while each surrounding sub-cell (i.e. Hi, ... , Hq) consists 
of one RS at the centre. For convenience, we approximate the coverage of BS and RS by hexagons as 
shown in Fig[TJ We define the central sub-cell (Hq) as base region and the six surrounding sub-cells 
(Hi, . . . , Hq) as relay region in every cell. We define a MS (user) present in the reference cell as 
target MS. We assume that BS and RS have Line of Sight (LoS) connection. All RSs are assumed 
to be decode and forward type relays, which decode the data received from BS and then forward it 
to the target MS. 

We consider universal frequency reuse, i.e. all cells use the same spectrum, which is shared 
between BS and six RSs. We consider the interference from the first tier of neighboring cells only. 
We also consider the effect of path loss and lognormal shadowing on the transmitted signal. Let 
BS transmits at power P to a MS located at distance d, then the received power at MS will be 
P(i~^10^/^°, where /3 is path loss exponent and ^ represents lognormal shadowing on BS-MS link. 
^ is a Gaussian random variable with mean dB and standard deviation a. Since thermal noise is 
negligible in an interference-limited reuse-one network, we ignore it in our computations. 

In practice, the association of MSs with BS or RS is determined based on SIR. If SIR experienced 
by MS from BS is above threshold, then it will be associated with BS otherwise with RS. In the 
present paper however, we have considered a simplified model where users present in BS region are 
associated with BS directly and those present in relay region are associated with the corresponding 
RS. It is assumed that the users are uniformly distributed in the respective regions of the cell. As 
explained later in Section |Vl out of the total call arrivals to the cell, a fraction is assumed to occur 

October 16, 2012 DRAFT 



4 



in BS region, while the remaining are assumed to have occured in relay region. 

We assume that K number of subcarriers are available in the reference cell, which are shared 
between BS and six RSs. Each RS and BS are allocated Krs and Kbs subcarriers respectively, such 
that K = Kbs + ^Krs. In this paper, we do not give any consideration to determine an optimal 
method for sharing the subcarriers amongst BS and six RSs. In the reference cell, all users have 
been allocated orthogonal subcarriers and therefore no intra-cell interference exists. However, in a 
network with universal frequency reuse, users will experience interference from RSs and BSs of 
neighboring cells. We analyze this interference on BS-MS, BS-RS and RS-MS link and compute 
their Cumulative Distribution Functions (CDFs). Using these CDFs, we determine the probability 
distribution of subcarrier requirement on these three links. We consider the rate requirement to be 
same for all users. The blocking probability on each link is calculated and then overall blocking 
probability for relay based OFDMA network is determined. 



In this section, we consider a target MS (user) in the reference cell. We analyze the SIR experienced 
by that target MS on BS-MS link if it is associated with BS and on BS-RS and RS-MS links, if it 
is associated with RS. Then, we compute the CDF of SIR on these links as derived in [?]. For this, 
we divide the incoming users into two groups, 

Group 1 : Users present in base region are associated with the BS directly on the BS-MS link. These 
calls are called direct calls. 

Group 2: Users present in the relay regions are associated with the BS via corresponding RS on 

BS-RS and RS-MS links. These calls are called hopped calls. 

Note that we use the terms users and calls interchangeably in this paper. 

A. SIR on BS-MS, BS-RS and RS-MS transmission links 

Let 'yss-Ms^ 1bs-rs and 'Jrs-ms denote SIR on a subcarrier used on BS-MS, BS-RS and RS-MS 
links respectively. Then we have. 



III. Inter-Cell Interference modeling 



PsMd^S-MS^^ 




IBS- MS — 



^iBS-MS ' 
10 



(1) 



Z^i=l ^BMdiBS-MS 



10 
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iBS-RS 



2^i=l ^BRO-iBS-RS 



PBRdss-Rs'^^ 




(2) 



and 



iRS-MS 



i=l ^RMO-iRS-MS 



PRMdj^S-MS^^ 




(3) 



where. 



• Pbm, Pbr and Prm denote the power transmitted by BS to target MS, BS to RS and RS to 
target MS respectively. 

• dss-Ms^ dss-Rs and d^s-Ms denote the distance between BS and target MS present in base 
region, BS and RS (with which the target MS is associated) and RS and target MS present in 
any of the relay regions respectively. 

• diBs-MS^ diBs-RS and di^s-MS denote the distance between i*^ neighboring BS and target MS 
(present in base region), i^^ neighboring BS and RS (with which the target MS is associated) 
and i*'* neighboring RS and target MS (present in any of the relay regions) respectively. 

• is the number of interferers in the first tier of cells. 

• ^bs-ms, iBS-RS and ^rs-ms represent lognormal shadowing on BS-MS, BS-RS and RS-MS 
links. Each of them is a Gaussian random variable with mean dB and standard deviations 
ctbs-ms, ctbs-rs and (Trs-ms dB respectively. 

• CiBS-MS, CiBS-RS and 6i?5-MS represent lognormal shadowing on i*^ neighboring BS and target 
MS link, i*'' neighboring BS and RS link, and i^'^ neighboring RS and target MS link. Each of 
them is a Gaussian random variable with mean dB and standard deviations OiBs-MS, cfibs-rs 
and (JiRs-MS dB respectively. 

B. CDF of SIR 

In this section, we determine mean and variance of interference to signal ratio Ibs-ms in two 
steps. In [?], [?] and [?], it has been argued that the total interference power received from various 

interferers (in a universal frequency reuse system) can be modeled by lognormal distribution with 
some mean and variance. We make the same assumption here. Accordingly, we proceed to calculate 
the mean and variance of Ibs-ms- 
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We rewrite Eq. [T] as, 



1 



1 



iBS-MS 



sr^N . d.Bs-MS N-^^Q ^'^g-"g-^^s-^^s IbS-Ms' 

Z-^t = l\ dgS-MS ' 



Step-1: Let (0, 0), (x, y) and (xj, Ui) be the coordinates of BS in reference cell, target MS in reference 
cell and the i^^ interfering BS present in the first tier respectively. 
Ibs-ms is grouped into two components, i?jS and C^s as. 



N 



IBS-MS 



BiCi 



(4) 



i=l 



where, Bi 



-/3 



2 , „ ^ „ SiBS-MS-^BS-MS 

and C, = 10 lo 



dBS-MS 

The position (x,?/) of target MS in the reference cell is random but the interfering BSs have fixed 
positions. Therefore, diBs-MS^ are correlated and as a result, BiS are correlated random variables. 

Ci is a ratio of two lognormal RVs, shadowing from i^'^ interfering BS to the target MS and 
shadowing from serving BS to target MS. As suggested in [?], it can be approximated by a lognormal 
RV with mean and variance + a^s-Ms)- Thus, all CjS are lognormal RVs but correlated. 

We assume that, 

E [d] = E [C,j\ and E [dCj] = constant, \f i^j. 

The first and second moments of Ibs-ms are determined as. 



E[Ibs- 



MS\ 



E [d] E 



N 



i=l 



(5) 



E [Ils-Ms] = E [C^] E 



N 



i=l 



+ E [ac, 



N 



N 



^5 Eb. 



,i=l 



J=l 



In Eq. [5] and [61 computations of E[Ci], E[Cf] and E J2i=i Bi are straightforward. E J2i=i E 



(6) 



1 2 



and E 



are solved as follows. 



Since the distances diss-MS^ between the target MS and the i^^ interfering fixed BS are correlated. 



E 



can not be separated into a sum of terms. It is computed by averaging over the area 
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E 



N 



-I 2 



3^3 



x,yeHo . 



N 

E 

i=l 



(x - Xj f + {y-yi 



dxdy. 



(7) 



Now, to compute E 



Sill ' expectation is taken over all possible positions (x, y) that target MS 



can take in the base region. These integrals are evaluated separately for each interfering BS and then 



summed for all BSs to get E J^Zi 



as shown below. 



E 



N 



1=1 



N 



N 



x,y€Ho 



{x - Xj)^ + {y- yif 
x"^ + y"^ 



-13 



dxdy. 



(8) 



Eq. U\ and [8] can be solved numerically for hexagonal geometry. 

Step-2: We have obtained first and second moments of Ibs-ms in Step-1 (Eq. [5] and [61). Its distribution 
can be approximated by lognormal distribution with parameters {i^Ibs-msi^^bs-ms^- 
In general the A;*^ moment can be written as. 



(9) 



Using k = 1 and 2 and on inverting, we obtain. 



1. 



2lnE[lBs-Ms\ + i^lnE[lBs-Ms\ 



and 



^Ls-Ms = [Ibs-ms] + InE [IIs-ms] 



Using Eq. [TO] and [IH we determine the distribution as, 

InX — fJ-lBS-MS 



(x) = $ 



X > 0. 



(10) 



(11) 



(12) 



^Ibs -MS 

Here $(x) is the standard normal CDF. 

Similar calculations are performed to obtain the CDF of Ibs-rs and Irs-ms on BS-RS and RS-MS 
links as, 

Unx - ,x>0 (13) 
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Inx — 



,x>0. (14) 



Thus, we have determined the distribution of interference to signal ratio on a subcarrier on the three 
transmission links, i.e., BS-MS, BS-RS and RS-MS links. 

IV. Analytical Model to determine resource requirement based on CDF of SIR 

In cellular OFDMA system, an incoming user is allocated a certain number of sub-carriers to 
satisfy its rate requirement. In our formulation, we consider that all incoming users have the same 
rate requirement R. Due to different SIR experienced by the users, they will require different number 
of subcarriers. 

The objective of BS is to satisfy the rate requirement of each user, by allocating it the requested 
number of subcarriers which depends upon its experienced SIR. There are Kbs orthogonal subcarriers 
available at the BS, each of bandwidth W Hz. 

Let, 'Jbs-ms experienced by a user while using m^^ subcarrier on BS — MS link. 

Then, the rate R achieved using M number of subcarriers on BS-MS link is given by, 

M 

R = Wj2iog,{l + jSs-Ms)- (15) 

m=l 

Since SIR on each subcarrier is assumed to be same, i.e. 7bs-m5 = 1bs-ms ■ ■ ■ Ibs- -MS ~ 1BS-MS^ 
the number of subcarriers (Af) required by any user can be expressed as, 

R-lo9io{2) ^^^^ 

W ■ logio {l + 'jBS-Ms)' 

Now, we distinguish users based on their sub-carrier requirement as follows. We divide the entire 
interference to signal ratio Ibs-ms range (determined in section UlI-BI) into / + 1 nonoverlapping con- 
secutive intervals with boundaries denoted by {iBS-MsYr^i- ^^^^ ^^"^ ^^^^ '^^ BS-MS link, when 
the received interference to signal ratio Ibs-ms f^Us in the range Ibs-ms ^ [^bs-ms^ ^bs^-ms]^ 
then user is considered to be in class r. As r lies in the range (1, ...,/ + 1) the highest possible 
class of a user will be / when Ibs-ms falls in the range Ibs-ms ^ [^bs-ms^ ^''bs-ms] ■ 

Let denote the number of subcarriers required by the class r user on that link. In the present 



October 16, 2012 



DRAFT 



9 

case, we have = r. Let A = ^logio (2). Then, Eq. [16] can be re-written as, 

Ibs^ms = (10^ - 1)-^ (17) 

Let, F BS-MS (M*") denote the probability that an incoming user belongs to class r and requires 
number of subcarriers on BS-MS link to meet its rate requirement. It is determined as, 

FbS~MS {M'') = P [IbS-MS < ^ BS-MS < F^V-Ms] (1^) 
= -^Ibs-ms {^bV-Ms) " -^Ibs-ms {^BS-Ms) ■ 

where, -^Ibs ms {^bs-ms) CDF of interference to signal ratio (Eq.[T2l) on BS-MS link. Similar 
calculations are performed to determine the probability of subcarrier requirement on BS-RS and RS- 
MS links and they are denoted by ^bs-rs (M^) and Frs-ms (M'') respectively. 

V. Analysis of Blocking Probability 

For a relay based cellular OFDMA system, we have two types of incoming calls (as mentioned 
in section Ull)) : direct and hopped calls. Let, No and Nh be the number of classes of direct calls 
and hopped calls respectively. Let, n and h denote the class of direct and hopped calls where, 
n = 1, . . . , Nf) and h = 1, . . . , N^- We denote the subcarrier requirement of class n of direct 
calls and class h of hopped calls by M|^, M^^^ and M^j^^ on BS-MS, BS-RS and RS-MS links 
respectively such that, = M^^^ + M^^^^.^. Note that denotes the total subcarriers required 
by the class h user of the hopped call. If an incoming user in base region requires M]^ subcarriers, 
then it belongs to class n of direct call and if it requires subcarriers in any of the relay regions, 
then it belongs to class h of hopped call. Now onwards, we will denote the probability of subcarrier 
requirement for various calls by Fbs-ms {Ml), ¥bs-rs {M^bb) and Frs-ms {M^rm) on BS-MS, 
BS-RS and RS-MS links respectively. 

To admit a direct call, the required number of subcarriers should be available at the BS. However, 
to accomodate a hopped call, the required number of subcarriers should be available at BS as well as 
RS. Thus, a direct call implies the arrival of one call on BS-MS link and a hopped call implies arrival 
of one call each on BS-RS and RS-MS link. In this section, we determine the blocking probability 
of users belonging to direct and hopped calls. 

We assume that call arrivals in each cell are Poisson distributed with mean arrival rate A. Let, a 
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fraction of the total call arrivals, say / be served directly by BS, then the arrival rate of direct calls 
is Ad = /A and that of hopped calls is A^ = (1 — /)A. The service times of each class of direct and 
hopped calls are exponentially distributed with mean K From the assumption of uniform distribution 
of users, hopped calls are equally distributed across the six RSs in the cell. Thus, the arrival rate 
of hopped calls in the coverage area of each RS is \h/Q- Let, P^^ and P^^ be the blocking 
probability of direct and hopped calls respectively. Then, the overall call blocking probability is 

Fij = /Pij, + (1-/)PBH- 

A direct call is blocked if the required number of subcarriers are not available at the BS and a 

hopped call is blocked if the required number of subcarriers are not available at any of the two i.e. 

BS or RS. We define the state of the system to be 

S = (M^f/^, . . . , M^t/S, . . . , M^^US^,MI,UI„ . . . , MM, • • • , M^"U^"), (19) 

where f/]^ and are the number of users in n^^ class of direct calls and /i*^ class of hopped calls 
respectively. denotes the number of subcarriers required by the n*^ class user of direct call and 
M'^j denotes the number of subcarriers required by the h^'^ class user of hopped call. 

These system states can be modeled by discrete time No + Nh dimensional Markov chain. The 
state space is finite and meets the following constraints- 

Nd Nh Nh 

J^M-U-+J2m^brU'h + J2^hrmU^ < Kss + Kns, (20) 

n=l h=l h=l 

Nn Nh 

Y^MIUI + Y^^hbrUI < Kbs, (21) 
n=l h=l 

Nh 

J2m'hrmU'h < Krs, (22) 

h=l 

> and > 0. (23) 

The constraint in Eq. [20] give an upper bound on the number of subcarriers available for allocation 
on the three links. The total number of subcarriers available at BS and corresponding RS gives an 
upper bound on the number of subcarriers that can be used in the system. The number of subcarriers 
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available at the BS gives an upper bound on the number of subcarriers that can be used by direct 
calls on BS-MS link and hopped calls on BS-RS link (Eq. [2T1) . Similarly, the number of subcarriers 
available at RS gives an upper bound on the number of subcarriers that can be used by hopped calls 
on RS-MS link (Eq. [22l) . In the system there can be either no user or a finite non-negative number 
of users on each link (Eq. [23l) . 

Example- r. For illustration, let us consider only one class of each call say, class 2 of direct call 
(i.e.. Ml = 2) and class 3 of hopped call (i.e.,M^ = 3). Let Kbs = 10 and Krs = 6. The 
number of subcarriers for hopped call {Mfj = 3) is the sum of subcarriers required on BS-RS and 
RS-MS links. Note that BS-RS and RS-MS link may require either Mfj^j^ = 1 and Mfjj^j^j = 2 or 
^HBR = 2 and Mjj^^^.j = 1 depending on the SIR experienced on each link. Thus, there are two 
possible combinations of subcarrier requirement for a hopped call on BS-RS and RS-MS links i.e. 
(1,2) and (2, 1). Let, the probability of subcarrier requirement of hopped call be P^^ (Mfj). Then 
we have, 

Fh {Ml = 3) = Fbs-rs (1) ^BS^MS (2) + Fbs-rs (2) Frs^ms (1) • (24) 
The arrival rate of hopped call of class 3 (Af^) and direct call of class 2 (A|,) is given by, 

\]j = XhFh [Ml = 3) , 

(25) 

-^i) = '^D^ BS-MS {M^d — 2) • 

The states of the system are represented by two dimensional Markov chain in Fig. [2l Each state 
corresponds to the number of subcarrier requirement for direct calls and hopped calls. There are 
various combinations of different subcarrier requirement on BS-RS and RS-MS links for a hopped 
call. 

The different combinations of subcarrier requirement for the first row of Markov chain in Fig. |2] 
denotes the case, when only users of hopped call are present. It is further illustrated in Fig. [3l where 
the state representation is modified to indicate the number of subcarrier requirement for hopped calls 
on BS-RS and RS-MS links distinctly. 

At any instant of time, the number of calls present in the system using various combination 
of subcarriers can be found by traversing a path as shown with dotted lines in Fig. [3l Similar 
combinations of subcarrier requirement of hopped calls with direct calls present in the system can 
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be obtained for various rows of Markov chain of Fig. [21 

As mentioned in the Example- 1, when a hopped call with Mjj = 3 arrives in the system, it 
requires either of the combinations (1,2) or (2,1) subcarriers on BS-RS and RS-MS links. This 
hopped call is blocked when the required number of subcarriers are unavailable at either BS or RS. 
Observing the dotted lines in Fig. [3l it becomes clear that after allocating the resources to 5^'^ user, 
BS is left with 1 subcarrier for new allocation on BS-RS link and RS has no subcarriers left for 
further allocation on RS-MS link. Similarly, other combinations of states leading to blocking state 
are shown in Fig. |3l This implies that there can be at most 5 users of hopped calls in this example. 

As can be noticed from this example, determining the set of all possible states which satisfy the 
given constraints for a single class of each call is complex. As the number of classes and the number 
of subcarriers at the BS and RS increase, the size of the state space increases and it becomes very 
difficult to determine all possible combinations. This complexity is due to two reasons: a) State space 
consists of the subcarrier requirement of calls of all classes on all the three links and b) The states 
on BS-MS and BS-RS links are interrelated because out of the Kbs resources available with BS, 
allocation for both direct and hopped calls is required to be done. 

To simplify the computational complexity, we consider the calls served by BS and RS as separate 
systems, because both have distinct set of subcarriers. For simplicity, we assume that for a hopped 
call, subcarriers are allocated independently by BS and RS. This is an important simplification which 
enables decoupling of state space of BS and RS. With this consideration, we determine the blocking 
probability in base and relay regions separately in the following subsections. The sum of blocking 
probability of calls in base region and relay region will give an upper bound on the overall blocking 
probability of the system. 

A. Blocking for Users present in Relay Region (Hopped Calls) 

When a user is in any of the relay regions and experiences SIR 'Jbs-rs and 'Jrs-ms on BS-RS 
and RS-MS link, it requires M^^^ number of subcarriers with probability ^bs-rs{Mhbr} ^^^^ 
M^RM number of subcarriers with probability ¥rs-ms{Mbiim)- The availability of subcarriers on 
both the links i.e. BS-RS and RS-MS links is determined. If subcarriers are available on both the 
links, M^Q^ and M^j^j^j subcarriers are allocated independently by BS and RS. Otherwise, that 
incoming user is blocked. In other words, blocking occurs when either M^br number of subcarriers 
are unavailable on BS-RS link or M^j^j^^ number of subcarriers are unavailable on RS-MS link. 
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Let, ^Bhbr '^Bhbm '^he blocking probabilty of hopped call on BS-RS and RS-MS link. 
Then, the average blocking probability of hopped calls (^Bh) given as, 

Pb, = 1 - (1 - Pb«,J(1 - PwJ. (26) 

In this subsection, we determine Pb^^^j,^ and in the next subsection we will determine Pb^^^. 

There are classes of hopped calls on RS-MS link, each requiring M^j^j^j subcarriers. The arrival 

rate of each h^'^ class of these calls at RS is = \H^B.s-AisiM^j^;^,^). Let the service time for all 

classes of call be exponentially distributed with mean service time ^. Then, the offered load for h^'^ 

a'' 

class on RS-MS link is p/i = It is assumed that after completion of a call, the subcarriers are 
released independently on both the links. 
We define the state of serving RS as, 

^RS = {M]j^mU\MI^^U\ . . . , M^^,,[/\ . . . , M^g^.f/^-), (27) 

where is the number of users of hopped calls of h^^ class and M^^^^^ is the number of subcarriers 
required by this hopped call of h^^ class. Any class of hopped call is said to be blocked, when all 
subcarriers Krs are in use. Therefore, the states of the system is modeled by Nh dimensional 
Markov chain. The state space is finite and meet the following constraints- 

Nh 

Y.M^rmU'' < Kns, U'' > and 1 < h < Nh. (28) 

h=l 

Example-2: Let us consider Krs = 4 subcarriers and Nh = 2 classes of hopped users. Let the users 
require Mjjj^j^j = 1 subcarrier with probability Fjis-MsiMjjj^j^j) = 0.6 and Mjj^j^^ = 2 subcarriers 
with probability ^rs-ms{Mhrm) = Thus, the arrival rate of class-1 users is = 0.6A and 
that of class-2 is (A^) = 0.4A. The states of the system are denoted by (M^^^jf/^ M^^jy^f/^). The 
state transition diagram is shown in Fig. IH Under the assumption of statistical equilibrium, the state 
probabilities are obtained by solving the global balance equations for each state. 

Let us consider any four interconnected states in Fig.lH If the flow in clockwise direction equals the 
flow in the opposite direction, then the process is said to be reversible [?]. Let, p{Mjjj^j^jU^ , ^hrm^"^) = 
p{l, 2) be the state probability. Then, from Fig. HI we have. 
Clockwise: 

p{l, 2) ■ ■ p{l, 0)-Xj,- p{2, 0)-Xl- p{2, 2) ■ 2/i. 
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Anticlockwise: 

2) • ■ p(2, 2) ■ yU ■ p(2, 0) ■ 2yU ■ 0) ■ Xjj. If these two expressions are equal then the process 
is said to be reversible [?]. 

We can express any state probability, say p{M}jj^j^jU^ , ^hrm^'^) terms of p{0, 0) by choosing any 
path between the two states, p(0,0) and the state itself, i.e, p{M}jj^jy.jU^ , Mfjj^j^jU'^) (Kolmogorov's 
criteria). 

In Fig. HI p(2, 2) can be obtained by choosing the path: (0,0), (0,2), (1,2) and (2,2), and we 
obtain the following equation, 

P{2,2) = ^,(^y-^-p{0,0) (29) 
2! V ^ / 

= ^-^-Mo.o). 

2! 1! ' ^ 

Thus, there are two users of class 1 and one user of class 2 and this state probability has product 
form. Similarly, Kolmogorov's criteria is applicable to a system with classes, and the state 
probabilities in A'^^^^ dimensional system will have product form [?]. 

Let, ^Urs probability that the system is in state ^Irs- Since all the states are reversible, the 

solution is given in the product form as, 

j-rh 

Pn,, = ^^'='^- (30) 

Let, Vth be the set of states in which an incoming hopped call on RS-MS hnk of either class is 
blocked. It is represented as, 

Nh 

Qh ■■= {^h eQRs-Yl ^hrmU^ > - M^hrm)- (31) 

h=l 

In Fig. m the states in which an incoming user of class-1 will be blocked are {(4, 0), (2, 2) and 
(0,4)}. The sum of the probabilities of these states is equal to the blocking probability for class-1. 
Similarly, the states in which an incoming user of class-2 will be blocked are {(3,0), (1,2) and 
(0,4)}. The sum of the probabilities of these states is equal to the blocking probability for class-2. 
The above illustration makes it clear that the blocking probability for any class can be obtained by 
summing the probabilities of all those states in which an incoming user of that class will be blocked. 

In general, for a relay based cellular OFDMA system with Nh classes, blocking probability for 



October 16, 2012 



DRAFT 



15 



hopped call of /i*'* class on RS-MS link is given by, 

^kuM - E^^- (32) 



h 



il 



The average blocking probability for hopped calls on RS-MS link is given by, 

Nh 

^Bhum = E ^BHUM^RS-MsiM^HRu)- (33) 

h=l 

In the next sub-section, we determine the blocking probability of direct calls. 



B. Blocking for Users present in Base Region (Direct Calls) 

When a user is in base region and experiences SIR -Jbs-ms, it requires number of subcarriers 
with probability ^bs-ms{M^). The availability of subcarriers is determined at BS. If they are 
available, then subcarriers are allocated by the BS. However Kbs subcarriers are also shared 
by h*'^ class of hopped call on BS-RS link, each of which requires M^^j^ subcarriers. Thus, there 
are Nd and Nh classes of direct calls on BS-MS link and hopped calls on BS-RS link respectively. 
The arrival rate of n*'^ class of direct calls is — bs-ms{Md) ^^^^^ of hopped calls 

is \)jBR = \h^bs-rs{MUj,). 

We define the state of BS as, 

^BS = {MhUh, MlUl, . . . , M^-U^-,M]j^^Ui, Ml^^Ul, . . . , M^%^U^"), (34) 

where U''^ is the number of direct users of n*'' class and is the number of hopped users of /i*^ 
class. M'^ and M^^j^ denote the subcarrier requirement of n*'* and /i*'^ class of direct and hopped 
calls respectively. If the subcarrier requirement for any class of direct call and hopped call on BS-RS 
link is same, then for BS both the calls will belong to the same class, irrespective of whether it is 
a direct or a hopped call. Thus, the state of BS can be modified as, ^bs — (MbsUbs) where m 
= 1,. . . , max(A^£), Nff), denoting the class of users arriving at the BS. M'^g denotes the subcarrier 
requirement m*'* class of user and U'g'g denotes number of users of m*^ class arriving at the BS. 

It is possible that some hopped calls get the required number of subcarriers on BS-RS link but 
not on RS-MS link. This is accounted by multiplying Xh by a discount factor 1 — ^Bhrm- Let, 
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the arrival rate of all calls at BS be Xbs- Then, the arrival rate of class m call at BS will be, 
^BS — + (1 ~ ^HRM^^H- The service time for all classes of calls at BS is exponentially 
distributed with mean service time ^. Then, the offered load for class m call at BS is pm = -jf-- 

Any class of calls (direct or hopped calls) at BS is said to be blocked, when all subcarriers Kbs 
are in use. Therefore, the states of the system are represented by max{Nu, Nh) dimensional Markov 
chain. The state space is finite and the constraints to be met are, 

max{N £, ,N u) 

^bsUbs < Kbs, U'Ss > ^ , I <m <max{ND,NH). (35) 

m=l 

Let, P^Bs be the probability that the system is in state VLbs- Since all the states are reversible, 
the solution is given in the product form as following [?], 

jjm 

Yimax{Nu,NH) {pmrBS 
^"ss ~ fAi M ^ f • ^^^^ 

ET-rmax{ND,NH) (pm) BS 
QBsLLm=l U^g\ 

Let, be the set of those states in which an incoming direct call or hopped call at BS of any 
class is blocked. It is represented as, 

max{N£),N[j) 

nm-^{^me^BS-- Yl MbsU'Ss>KbS-M^s}- (37) 

m=l 

Therefore, the blocking probability for m*'* class user at BS is given by, 

PS,. = E^Ob. (38) 



= E 



nmax{ND,NH) {Pm) BS 
m=l U^g\ 



ymax{ND,NH) [pm) BS 



Qbs llm=l Vggl 

The average blocking probability for direct calls on BS-MS link (Ps^,) and hopped calls on BS-RS 
link {^Bhbr) is given by. 

No 

^Bo-Y.^Bbs^bs-ms{M^), 

(39) 



^BBBn -Y.^'^bs'^BS-Rs{M^Br)- 



m=l 
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VI. Numerical Results 

In this section, we present results based on the analytical models developed in the previous sections. 
We consider the downlink of relay assisted OFDMA system. The values of system parameters chosen 
for the analysis are as per the LTE standard [?] and are given in Table HI We do the analysis for 
four values of the rate required by the users - 64, 128, 256 and 1024 Kbps. 

Fig. |5] and Fig. |6] give the probability of subcarrier requirement for four different data rates on 
BS-MS and BS-RS link respectively. It is the probability of a call belonging to a certain class. 
The probability of a call belonging to lower class is more on BS-RS link due to line of sight path 
and lesser impact of shadowing. From these two figures, we observe that the subcarrier requirement 
changes with the change in rate requirement (R). For lower R, less number of subcarriers are required 
and therefore, majority of users will belong to lower classes. For example, in Fig. [51 for i? = 64 
Kbps, the probability of a call belonging to a certain class is non-zero from class 1 to class 15 and the 
probability of call belonging to higher classes is close to zero. Similarly, for high R, say R = 1024 
kbps, there are effectively no users that require less number of subcarriers. Therefore, the probability 
of a call belonging to a certain class is non-zero from class 5 to class 35 and the probability of 
call belonging to lower classes is almost zero. This clearly indicates two things- that for a given 
rate requirement, 1) the number of class will depend on the range for which the probability of a 
call belonging to a certain class (i.e., subcarrier requirement) is non-zero and 2) the definition of 
class will not always be iVf = r (Section HVT). On the contrary, it will depend on the lowest and the 
highest class for which the subcarrier requirement is non-zero. For example, for R = 1024 Kbps, 
there will be 31 classes and it will range from class 5 to class 35, i.e., = r + 4 (offset of 4). 
Thus, definition of class and determining the number of classes depends on the rate requirement. 

Fig. |7] and Fig. [8] illustrate the impact of rate requirement on blocking probability for two cases: 
when the subcarrier bandwidth considered is of 15 KHz. and 30 KHz. We observe that increasing 
the subcarrier bandwidth results in an increase in the blocking probability. It is because when the 
subcarrier bandwidth is more, the total number of subcarriers available in the system reduces. In 
this case, even though a user may meet its rate requirement with fewer number of subcarriers, the 
overall blocking probability is likely to increase. Similar observation can be made from Fig. |9] where 
the blocking probability is computed for three different subcarrier bandwidth, 15, 30 and 60 KHz, 
for fixed rate requirement of 64 Kbps. 
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From Fig. |7] and Fig. [8l we also observe that irrespective of the subcarrier bandwidth, the blocking 
probability is influenced by the rate requirement of users. As the rate requirement of users increase, 
they will require more number of subcarriers and therefore, blocking probability increases. 

As an intuitive insight, when the subcarrier bandwidth is high, the blocking probability is influenced 
by the number of users belonging to lower classes and higher classes. If majority of users belong to 
lower classes, then the subcarrier bandwidth will result in an allocation which will be much more 
than their requirement, leading to an inefficient resource utilization and hence, increase in blocking 
probability. On the other hand, if majority of users belong to higher classes, then they will quench 
their resource requirement in fewer resources and blocking probability is likely to reduce. However, 
in general when the probability that a user belongs to a higher or a lower class (i.e., the subcarrier 
requirement of a class is more or less) is equal, an increase in the subcarrier bandwidth will reduce 
the number of resources available in the system and hence there will be an increase in the blocking 
probability. 

VII. Conclusions and Future Work 

In cellular OFDMA networks, in order to meet the same rate requirement, the number of subcarriers 
required are different for different users (due to differences in their locations and experienced SIR) 
on various links. Therefore traditional methods of blocking probability computation cannot be used 
directly. We have proposed an analytical model to evaluate blocking probability for relay based 
cellular OFDMA networks. The CDF of SIR is determined and is used to compute the probability 
distribution of subcarriers required on the three transmission links. The incoming users are classified 
into different classes based on their subcarrier requirement. We have modelled such a system by a 
multi-dimensional Markov chain. The effects of subcarrier bandwidth (W) and rate requirement (R) 
on the blocking probability are analyzed. We have also analyzed the effect of rate requirement on the 
definition of class and number of classes. This model is useful in the performance evaluation, design 
and planning of relay based cellular OFDMA networks like LTE. In this work, we have assumed 
that the total subcarriers are partitioned between BS and RS in somewhat heuristic manner. We 
are investigating the impact of optimal resource partitioning between BS and RSs on the blocking 
probability as future work. 
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TABLE I 

Parameters for Numerical Analysis 



Parameters for Numerical Analysis 


Value 


Inter BS distance (meters) 


1732 


System Bandwidth (MHz) 


10 


Subcarrier bandwidth (W) (KHz) 


15 (30) 


Number of Subcarriers available at 


480 (when 


BS (Kbs) 


W = 1^ KHz) 
240 (when 
= 30 KHz) 


Number of RSs in each reference cell 


6 


Number of subcarriers available at 


30 (when 


each RS (Krs) 


W = lh KHz) 

15 (when 

1^ = 30 KHz) 


Number of mtertemng cells from first 




tier (iV) 


6 


Path loss exponent 


3.5 


Shadowing standard deviation on 




BS-MS link {ass-MS = cTiBs-Ms) 
and RS-MS link {aB^s-MS = (^ws-ms) 


8 dB 


Shadowing standard deviation on 




BS-RS link {ass-RS = cFiBS-rs) 


4 dB 


Rate requirement of each call {R) (Kbps) 


64, 128, 
256, 1024 


Maximum number of classes for each call 




{Nd = Nh) 


10 (15) 


Fraction of calls arriving at BS (/) 


0.5 


Mean Arrival Rate (A) (calls/unit time) 


1 to 80 
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Fig. 1. Architecture of Relay based Cellular OFDMA System 
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Hopped Call 




Fig. 2. Two-dimensional state transition diagram for a system with Kbs = 10, Kjis = 6, = 2 and Mjj = 3. Each call has 
only one class. 



October 16, 2012 



DRAFT 




October 16, 2012 



DRAFT 



23 



Class 1 of hopped call 




Fig. 4. Two- dimensional state transition diagram for a system with Krs = 4 and Nh = 2. Incoming users are divided into 2 
classes 




Fig. 5. Probability of Subcarrier requirement (i.e. class) for varying rate requirements on BS-MS link (Subcarrier BW = 15 KHz) 
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5 10 15 20 25 30 35 
Class Number 



Fig. 6. Probability of Subcarrier requirement (i.e. class) for varying rate requirements on BS-RS link (Subcarrier BW = 15 KHz) 




Mean Arrival Rates 



Fig. 7. Impact of Rate Requirement on Blocking Probability when Subcarrier BW = 15 KHz (Number of classes for 64 and 128 
Kbps = 15, Number of classes for 512 and 1024 Kbps = 31) 
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Fig. 8. Impact of Rate Requirement on Blocking Probability when Subcarrier BW = 30 KHz (Number of classes for 64 and 128 
Kbps = 15, Number of classes for 512 and 1024 Kbps = 31) 




Fig. 9. Variation in Blocking Probability with change in Subcarrier BW for a given rate requirement = 64 Kbps and Number of 
classes = 15 
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